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Phase diagrams of the one-dimensional extended Hubbard model (including nearest-neighbor inter- 
action V) at half- and quarter-filling are studied by observing level crossings of excitation spectra 
using the exact diagonalization. This method is based on the Tomonaga-Luttinger liquid theory 
including logarithmic corrections which stem from the renormalization of the Umklapp- and the 
backward-scattering effects. Using this approach, the phase boundaries are determined with high 
accuracy, and then the structure of the phase diagram is clarified. At half-filling, the phase diagram 
consists of two Berezinskii-Kosterlitz-Thouless (BKT) transition lines and one Gaussian transition 
line in the charge sector, and one spin-gap transition line. This structure reflects the U(l) ® SU(2) 
symmetry of the electron system. Near the U = 2V line, the Gaussian and the spin-gap transi- 
tions take place independently from the weak- to the intermediate-coupling region, but these two 
transition lines are coupled in the strong-coupling region. This result demonstrates existence of a 
tricritical point and a bond-charge-density-wave (BCDW) phase between charge- and spin-density- 
wave (CDW, SDW) phases. To clarify this mechanism of the transition, we also investigate effect of a 
correlated hopping term which plays a role to enlarge BCDW and bond-spin-density-wave (BSD W) 
phases. At quarter-filling, a similar crossover phenomenon also takes place in the large-!/ region 
involving spin-gap and BKT-type metal-insulator transitions. 
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I. INTRODUCTION 

One-dimensional (ID) electron systems have been ex- 
tensively studied motivated not only by theoretical in- 
terest but also by the discovery of quasi-lD conductors 
and high-T c superconductivity. In the ID electron sys- 
tems, due to the charge-spin separation, the low-energy 
excitations in the charge and the spin sectors may have 
gaps independently, and then various phases can appear. 
However, phenomena caused by interplay between these 
two degrees of freedom have not been fully understood 
even in simple models. In this paper, we turn our at- 
tention to the phase transitions in the so-called extended 
Hubbard model (EHM), 

Hehm = -t^2(c\ s c i+liS +H.c.) 

is 

+U y] n^nn + V riim+i, (1) 

i i 

at half- and quarter-filling, where both charge and spin 
gaps can open. 

The EHM at half-filling has been studied using var- 
ious approaches. In the weak-coupling limit, J.hxjL.ahase 
diagram is analytically obtained by the g-ologyllotra (see 
Appendix A). According to the result, there appear in- 
sulating charge- (CDW) and spin-density-wave (SDW) 
phases, and metallic phases where the singlet supercon- 
ducting (SS) or the triplet superconducting (TS) corre- 
lation is dominant [see Fig. l](b)]. On the other hand, in 
the strong-coupling limit, the perturbation theory gives 



the phase boundary of tbjajQDW-SDW transitionElErO and 
of the phase separationJj'LrllJ Thejpsj^f-jegkin has been 
discussed by numerical analysis. BEjOooIIj However, 
the phase diagrams are not fully understood, because 
the charge and the spin gaps open exponentially slow [see 
Eqs. (18) and (pl|)], which makes it difficult to determine 
the phase boundaries by the conventional finite-size scal- 
ing method. Especially for the transition between the 
CDW and the SDW phases, even the property of the 
transition itself is not clear, because the transition is of 
the second order in the weak-coupling theory, while it is 
of the first order in the strong-coupling theory. 




(a) X < (b) X = (c) X > 

FIG. 1. Weak-coupling phase diagrams of the EHM includ- 
ing the correlated hopping term (Ref. |o|). The phase dia- 
grams are given by combinations of a Y-shaped structure for 
the charge part [two BKT lines (solid) and one Gaussian line 
(dashed)] and an I-shaped one for the spin part. 

Recently, the author has. clarified the mechanism of 
the CDW-SDW transition£3 According to the result, the 
phase boundary consists of two independent transition 
lines, and the crossover of the CDW-SDW transition is 
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related with whether these two transition lines are sepa- 
rated or coupled. The result also demonstrates the exis- 
tence of the bond-charge-density-wave (BCDW)E-O phase 
in the very narrow region between the CDW and the 
SDW states. In this paper, we not only give the details 
of the letter, but also clarify the entire phase diagram of 
the EHM at half-filling. 

In order to clarify the above scenario for the phase 
transition between the CDW and the SDW phases, we 
also consider generalizing the EHM b^addias the fol- 
lowing correlated hopping interactions:EJ'l 19 lr°rEj 

Hx = X ^(ct s c i+1 , s + H.c.)K_ s - n i+ i,_ s ) 2 . (2) 



This interaction can be derived ras a site-off-diagonal el- 
ement of the Coulomb integralJla Especially, the three- 
body part is justified as an effective interaction of the 
three-band model.EJ The weak-coupling phase diagram 
is known by the g-ologycj as shown in Figs. [I] (a) and (c) 
(see Appendix A). In these phase diagrams, the two tran- 
sition lines between the CDW and the SDW phases do 
not synchronize. And a BCDW or a bond-spin-density- 
waveE3 (BSDW) phase appears. The analysis of the gen- 
eralized model will clarify the tricritical behavior in the 
pure EHM. The final results are shown in Fig. [l^. 

A charge-gap phase is known to appea r n n o tj only 
at half-filling but also at quarter-filling,E3'tSEaj2j due 
to the-jcfflhetpof the Umklapp scattering in the higher 
order .Ej EaSEl In this case, the interplay between charge 
and spin instabilities is also expected. In fact, we will 
conclude that a crossover phenomenon also exists in the 
large- V region at quarter-filling (see Fig. |l6| ). 

Throughout this 
paper, we use the le\jfeL- [Cr . pssin | g . sroproa ch to determine 
the phase boundaries.HME3SaOEaO This method is. 
based on the Tomonaga-Luttinger (TL) liquid theoryo 
(which is equivalent to the c = 1 conformal field theory) 
including the logarithmic corrections, which stem from 
the renormalization of the Umklapp- and the backward- 
scattering effects. In the theoretical scheme, the tran- 
sition points are identified by the level crossing of the 
excitation spectra in the finite-size ring with size depen- 
dence 0(L~ 2 ), where L is the system size. Therefore, the 
phase boundaries are obtained with high accuracy, using 
the numerical data of finite-size clusters. 

This paper is organized as follows. In Sec. [□], we re- 
view the level-crossing approach based on the TL liq- 
uid theory and the renormalization group, developed in 
Refs. |3H-|3r 



In Sec. Ill 



we discuss the discrete sym- 
metries of wave functions to connect excitation spectra 
and the corresponding physical states. In Sec. IV 



we 



discuss the character of the phases that appear in the 
phase diagrams. In Sec. [v], we analyze the instabilities of 
the EHM at half-filling, and clarify the phase diagram. 
In Sec. VI, we analyze the metal-insulator transition at 
quarter-filling. Finally, a summary and discussions are 
given in Sec. VII. In Appendix A, we briefly explain the 



traditional g-ology analysis for the generalized EHM at 
half-filling. 



II. PHASE BOUNDARIES 

First, let us perform a general argument fox 
tron systems based on the bosonization theory.Er 
The continuous fermion fields are defined by Cj 

4>L,s(x) 

with 



■0R,s(^) (the lattice constant a — > 0, x = ja) 



1pr,s( x ) 



U, 



Tjs irk F x i/V2-\r(^ +s<t,„)-e -s9„] 



(3) 



where r = R, L and s =t> I refer to -I- and — in that 
order, a is a short-distance cutoff, kp is the Fermi wave 
number defined by kp = 7rn/2, with n being the electron 
density. The field <j> v and the dual field Q v of the charge 
(y = p) and the spin (y = a) degrees of freedom satisfy 
the relation 



[<^(x), 0u(x')] = -ytf^signfx - x'). (4) 

U r ^ s ensures awtieommutation relations of the different 
fermion fields. r 8 ^ 40 ! These operators are Hermitian and 
satisfy the relation 



(5) 



Using the formalism, the low-energy behavior of ID elec- 
tron system can be described by the sine-Gordon models 
for the charge and the spin sectors. When 2q (q: integer) 
electrons contribute to the Umklapp scattering, the ef- 
fective Hamiltonian for the system with length L is given 
by 

n = E £ f &x[K„{d x e„) 2 + k- 1 ^) 2 } 



2gi± 
(27ra) 2 J, 

2fl 3 ± 
(27ra) 2 

2 ff3|| 
(2na) 2 J 



dxcos[v / 8^ IT (a;)] 

dx cos[qV8<j) p (x) + 5x] 

dxcos[qV8(j)p(x) + Sx] cos[v / 8'/' -(a;)], (6) 



where v v and K v are the velocity and the Gaussian cou- 
pling, respectively for each sector. and denote 
the amplitude of the backward and the Umklapp scat- 
tering, respectively. The Umklapp term vanishes except 
for the case 5 = 2pir — 4qkp = where p is also an 
integer, and p/q is an irreducible fraction. Thus, the 
electron filling that a charge gap-caa-paen is quantized 
to commensurate cases n — p/qtJ^BcB This condition 
can also be deru£ed-,from the generalized Lieb-Schultz- 
Mattis theoremE3'E2l In this paper, we consider q = 1 
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(half-filling) and q = 2 (quarter-filling) cases. At half- 
filling, in the weak-coupling limit, the couplings of the 
backward and the Umklapp scattering for the EHM are 
identified as gi± = —gs± = U — 2V (see Appendix A). 

In addition to the g%± term, there exists another Umk- 
lapp operaioiiwith coupling constant <? 3 || , which transfers 
finite spin.olI3'E3 In the weak-coupling limit, this param- 
eter is identified as <7 3 || = —2V. In the present analysis, 
we will not consider this term explicitly, because the scal- 
ing dimension of this term is always higher than that of 
the other nonlinear terms in Eq. (^J). However, in the 
strong-coupling region, a charge-spin coupling effect may 
appear due to this term. 

If these nonlinear terms are absent (gi±_ = g-3±_ = <?3|| = 
0), the excitation spectra and their wave numbers in the 
finite-size system are described by 



where O u o is the "marginal field" ,0 and the deriva- 
tives are defined by d, B = (v^df =F \d x )/2 with imag- 
inary time t. This operator corresponds to particle- 
hole excitations near the right and the left Fermi points 
(N u = N v = I). O v \ and O v i are linear combinations of 
current excitations {m u — ±q). O v z is an excitation ac- 
companying variation of the number of electrons or spins 
(n„ = ±1). We have to choose antiperiodic boundary 
conditions ipr,s(x + L) — ~ip rtS (x) to extract the exci- 
tation spectra for.. Q v \ and V 2 fields, when q is odd, 
and 0„3 field,oE3c3 because ths phase fields satisfy the 
following boundary conditions :E3 



b v {x + L) = 4> y (x) ~ V2im v , 
9 v (x + L) = v {x) + V27rm„, 



(14a) 
(14b) 



E — Eq — 



2ttv p 
2?r 



2ttv„ 



P - P = —(s p + s a ) + 2m p k F , 



(7) 
(8) 



where the scaling dimensions and the conformal spins are 
given by 



and the Fermi operator is given by these phase fields as 
in Eq. (|). 

The effects of the gi± and the g%± terms in Eq. (g) are 
renormalized in the scaling dimensions logarithmic 
corrections which are analyzed by the renormalization 
group (RG) equations derived under the change of the 
cutoff a — ► e d/ a (Ref. ^2j). Within the one-loop order, 
the RG equations are given by 



Xv = \ (jit + m » K ^) +N » + N *i 

s v = n v m v +N U — N v . 



(9) 
(10) 



Here n p is the change of 2n p electrons, and n a is the 



d^o(Q 

dl 
dy V( p{l) 

dl 



-yvo{l)yv<p(l), 



(15a) 
(15b) 



total z-spins = n a . 

particles moved from the left charge (spin) Fermi point 
to the right one. The non-negative integers N v and N v 
are the particle-hole excitations near the right and the 
left Fermi points, respectively. The scaling dimensions 
are related to the critical exponents for the correlation 
functions in the large distance as 



(Oi(r)Oi{r')) oc |r - r '\- 2 ^ i+x "^ 
where the operator is given by 



W denotes the^ number of where yp0 ( ) = 2{q 2 K p - 1), y a0 (0) = 2{K a - 1), 

y P 4>{ti) = 93J-/nv Pl y<j<t,(0) = gi^/irvcr, and we have set 
I = InL. These equations determine the RG flow dia- 
grams. Note that there is a difference between the cases 
for the charge and the spin sectors reflecting their sym- 
metries (see Fig. ||). In the following subsections, we 
discuss the phase transitions for each sector described by 
these RG flow diagrams. 



(11) 



n. = np n° 

O v n ^ mv EEe iv/5(n ^ +m '^" ) , (12) 

or linear combinations of these operators. Therefore, 
there are one-to-one correspondences between the exci- 
tation spectra and the operators. 

Now we turn our attention to the excitation spectra 
which correspond to the following operators: 

O v0 = -^Bfadcpv, (13a) 

O vX = V2cos(qV2<t> u ) oc O u aq + (D» _ q , (13b) 

O v2 = V2~sm(qV2(f>„) oc 0^ q - 0^ q , (13c) 

O v3 = exp(iV26 v ) = Olo, (13d) 



A. Spin-gap transition 

First, we consider the phase transition in the spin de- 
gree of freedom (y = a). The spin sector with an SU(2) 
symmetry belongs to the universality class of the leveL-1 
SU(2) Wess-Zumino-Novikov-Witten (WZNW) modelS 
In this case, the RG flow in Fig. || (a) is fixed on the 
y<jo{l) = y<T4>{l) line. Then for y a o{l) > 0, the exponent 
is renormalized as K* = 1 , and the solution of Eq. jl5| ) 
is obtained as 



y<ro(l) 



y CT o(0)/ + l' 



(16) 



where y CT o(0) is the bare coupling constant. Combin- 
ing the renormalized coupling and the operator-product- 
expansion coefficients, the singlek-|k&jrf-} and the triplet 
(Xo-2,3) excitation spectra split asEJ'C3Eil 
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1 3 

£<xi(0 = 2 + J^oCO) 

^2,3(0 = 2 - |^o(0- 



(17a) 
(17b) 



When 2/o-o(0) < 0, yo-o(l) is renormalizcd as y,xo(^ ^ 
oo) = — oo, then a spin gap appears. At the critical 
point [yo-o(O) = 0], there are no logarithmic corrections in 
the excitation spectra. Therefore, the critical point is ob- 
tained by the intersection of the singlet and the triplet ex- 
citation spectra [x a \ — £ CT2 ,3)-"i25t3ej This level cross- 
ing corresponds to the condition for the spin-gap phase 
boundary gi± = g a = in the standard g-ology analysis, 
(see Appendix A and Table ||) . 




FIG. 2. RG flow diagram for (a) the spin and (b) the 
charge degrees of freedom. For the spin sector, the RG flow 
is fixed on the y a o = y<j4> line due to the SU(2) symmetry, 
and the spin-gap transition takes place at j/ CT o = Va^, = 0. 
For the charge sector, BKT-type transitions take place on 
the y p o — iy p< (, lines with y p o > 0. When q — 1, the 
j/po = +y P 4> line corresponds to the SU(2) symmetry of the 
jy-paring. The y p o = —y P <j, line reflects the hidden symmetry 
of the sine-Gordon model. A Gaussian transition occurs on 
the y pr f, = line with y p0 < 0. 

The asymptotic behavior of the spin gap against a pa- 
rameter of a model A near the critical point A c is obtained 



by the two-loop RG equation ani 
lation length y cr0 (ln^) ~ —1 ast 



the definition of corre- 



Acr ~ Vtr/i oc - A c exp[-const/ (A - A c )], (18) 

where we have used a relation A — A c oc 1 2/crO (0) | . Note 
that this is the same behavior as that of the spin-gap in 
the ID negative-C/ Hubbard model at half-fillingo 



B. Berezinskii-Kosterlitz-Thouless transition 

1. SU(2) symmetric case 

Next, we consider the instabilities in the charge sector 
(y = p), which are described by the RG flow diagram 
given in Fig. |(b). At half-filling (q = 1) and V = (the 
Hubbard model) , the sign of the on-site interaction U in 
the Hamiltonian_p) is reversed by the following canonical 
transformation:!!]! 



(-l)V 



(19) 



This transformation also projects the spin (('-pairing) op- 
erators onto the ^-pairing ones 

4 = (~l)%4v ^ = {- l T^lc lh rfi = \{n, - 1), 

(20) 

without losing the SU(2) svmmetryEl It follows from 
Eq. (|J) this transformation corresponds to the replace- 
ment of the indices as p <-> a. Therefore, the spin part 
of the sine-Gordon model of Eq. (0) is mapped onto the 
charge part, and the operators O p \ and P 2, P 3 denote 
the "singlet" and the "triplet" for the charge part, respec- 
tively. Thus, the exponent is renormalized as K* = 1 for 
U < 0, and the charge gap opens for U > 0. 

In the case when the SU(2) symmetry in the charge 
sector is broken by finite V, O p i, O fi 2, and P 3 refer 
to "dimer", "Neel", and "doublet", respectively, by fol- 
lowing Ref. [32]. Then, if the initial value of the RG 
flow moves across the SU(2)-symmetric line [y P< />(0) = 
y P o(0) > 0], a_Hep^zinskii-Kosterlitz-Thouless (BKT)- 
type transitiorO'E3c3 occurs between the TL liquid phase 
and the twofold-degenerate gapped state. For tiiis tran- 
sition, one can show that a charge gap opens asc3 



4 



Ap oc exp(— const/y^ A — A c ), 



(21) 



-* <t> P + tt/v8. 



This operation interchanges the roles of the operators 
O p i and P 2 as 



COS V^p 

sin V2(f) p 



— sin \i2<j> p , 
cos V2(f)p . 



(23a) 
(23b) 



Therefore, this symmetry indicates that the SU(2)- 
symmetric line in the RG flow diagram is mapped onto 
the opposite side of the y p o axis, and another BKT tran- 
sition may occur at y p( j,(0) — —y p o(0) < 0. We call the 
symmetry of this BKT line "hidden SU(2) symmetry" . 
Since this symmetry originates from that of the sine- 
Gordon model, it is not contained explicitly in the orig- 
inal Hamiltonian. The renormalized scaling dimensions 
of O p i , P 2 , and P 3 near the critical lineu^f-the hidden 
SU(2) symmetry are calculated as follows:c3t3 



Xpi(l) = l-\y P o(l)(i + 2t), 



x p2 {l) = i + 7^0(0(1+ |*) ; 



1 3 

2 ' 4 

WO = 7, ~ T2/po(0 



(24a) 
(24b) 
(24c) 



Therefore, the critical point for the hidden SU(2) BKT 
transition can be determined by the level crossing be- 
tween the "dimer" .and the "doublet" excitation spectra 



pi 



L p3 



< X 



P2 



).E3 This level crossing corresponds to 



< in the g-ology analysis (see 



the condition g 3 ± = —g 
Appendix A and Table 



2. Non-SU(2) symmetric case 

We consider the BKT transition for q > 2 case. This 
situation may appear in the metal-insulator transition at 
quarter-filling (q — 2). In this case, the critical line no 
longer has an SU(2) symmetry. However, by replacing 



the variables as 



9' = 9 p /q, and K' 



q 2 K n 



where |A — A c | oc t, and t = \y P 4>(l)\/ypo(l) — 1 stands 
for the deviation from the BKT critical line. Note that 
Eq. (|2l]) is a different asymptotic behavior from that of 
the spin-gap transition described by Eq. (|is|), so that 
we discriminate the spin-gap transition from BKT-type 
transitions in this paper. The critical point for this BKT 
transition can be obtained without calculations, because 
it is fixed by the SU(2) symmetry of the Hamiltonian. In 
the case of the EHM, the BKT transition line is fixed on 
the V = line for U < 0. 

Now, we consider the region for y p< p(l) < 0. The sine- 
Gordon model has a symmetry under the transforma- 
tion to reverse the sign of the nonlinear term cos\/80p- Xp^/q < x p i 
This transformation corresponds to the shift of the phase 
fields: 



pi 

the sine-Gordon model for the charge part of Eq. (|6|) is 
mapped onto the case of q = 1. Then, the BKT tran- 
sition between the TL liquid and the 2g-fold-degenerate 
gapped state takes place when the renormalized expo- 
1/q 2 . The scaling dimensions for 



nent becomes K* = 

the Op\ and the O a 2 fields near the BKT critical line 



'pi 



remain unchanged, while the O p s field changes as 



x p3 (l) 



(25) 



Therefore, the BKT critical point corresponding to y p $ 
y p o > is given by the level-crossing point of x P 2 



There is another excitation spectrum that can be 
used to determi ne th e BKT critical point. This is the 
"marginal field'j- 



(22) sion is given b> 



13a) whose renormalized scaling dimen- 



ajpo(0 = 2-tfpo(0(l + |*) 



(26) 



In this case, the critical point can be determined by the 
level crossing of x p o = 4cc P 2. 



C. Gaussian transition 

In addition to these BKT-type transitions, a Gaus- 
sian transition occurs at y p( p(0) — and y p o{0) < 0. 
This is a second-order transition between the two gapped 
states which corresponds to the different fixed points 
[Vp4>{1 — > oo ) = ioo]; an d the gap vanishes just on the 
critical point. The transition point is given by the level 
crossing between the "dimer" and the "Neel" excitations 
(xpi = x P 2 < x P 3), because the O p \ and the P 2 fields 
inte rchang e their roles at y p <p{0) — as was explained in 
Sec. II B 1 . In the g-ology, this level crossing corresponds 
to the condition — with g p < (see Appendix 
A and Table ||). Since the nonlinear term vanishes on 
the critical line, there is no effect of the renormalization. 
Therefore, the scaling dimensions on the Gaussian line 
are given by 



Xpi 
Xp3 



Xp2 - 
1 

2Ka 



(27a) 
(27b) 



without logarithmic corrections, and K p < 1 is satisfied. 
The asymptotic form of the gap near the Gau ssian tran- 
sition point can be obtained by solving Eq. (15b) with 
an approximation y p o{l) ~ y p o(0) |£und definition of the 
correlation length |y p<p (ln ^)| ~ 1 as£3 



v p /£, oc |A - A c | ! 



(28) 



where we have used the relation y p< p(0) oc A — A c . The 
valley of the gap becomes steeper as K p is decreased. 
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g-ology 


level 


crossing 


spin gap 


ffiJ-(= 9*) = 




= Xa2,3 


SU(2)BKT 


S3i = 5p > 


Xpl = 


X p3 ^ 3<pl 


hidden SU(2)BKT 


93± = -flp < 


Xpl = 


£p3 < a; P 2 


Gaussian 


53± = 0, 5 P < 


Xpl = 


X P 2 < X p 3 



TABLE I. Correspondence between the g-ology and the 
level-crossing approach at half-filling. The scaling dimen- 
sions Xvi correspond to excitation spectra of singlet (x a i), 
triplet(a; -2,3), "dimer"(a;pi), "Neel" (xp2), and "doublet" (x P 3 
states. Examples of these level crossings are shown in Figs 
and ^. 



III. DISCRETE SYMMETRIES 

To perform the level-crossing analysis discussed in 
Sec. O, we need to identify the relevant excitation spec- 
tra. For this purpose, we discuss the discrete symmetries 
of wave functions corresponding to the excited states. 
The physical meaning of these symmetries will be clari- 



fied in Sec. IV 



The discrete symmetries are defined under particle- 
hole (C : Ci S — > (— l) 4 c- s ), space-inversion (V : c JS — > 
CL-i+x,s)i an d spin-reversal (T : Ci S — > Ci— S ) transfor- 
mations. They give eigenvalues ±1. In addition, shift 
operation by one site (S : Ci S — > q+i, s ) is defined which 
has an eigenvalue e lk . The symmetries of wave functions 
can be explained by combining those of the ground^tate 
and those of the operators for the excited states.E3 For 
the ground state of the EHM, we choose periodic (an- 
tiperiodic) boundary conditions when N/2 is odd (even), 
where N is the number of electrons. Then, according to 
the Perron-Frobenious theorem, the discrete symmetries 
of the ground state are C = V = T = 1 and k = 0, if 
we choose the representation for the basis and use the 
symmetry operations defined in Ref. Ba 





operators 


C 


V 


T 


k 


BC 


G.S. 


1 


1 


1 


1 





±1 


marginal 




1 


1 


1 





±1 


CDW 


sin V2~4> p ■ cos V2(j>cr 


-1 


-1 


1 


2k F 


±1 


SDW Z 


cos V2<f>p ■ sin \/20ct 


-1 


-1 


-1 


2fc F 


±1 


SDW± 


cos y/2<j) p ■ exp ±i\/29 a 


* 


1 


* 


2k F 


±1 


BCDW 


cos V2~4> p ■ cos V2~4><t 


1 


1 


1 


2k F 


±1 


BSDW Z 


sin yf2(f>p ■ sin \f2<j) a 


1 


1 


-1 


2k F 


±1 


BSDW± 


sin \/2~4>p ■ exp ±\\/29 a 


* 


-1 


* 


2k F 


±1 


ss 


exp \\/2Qp ■ cos V2~4>a- 


* 


1 


1 





±1 


TSo 


exp \\/2Qp ■ sin \/2(\> r y 


* 


-1 


-1 





±1 


TS±i 


exp \\p2Qp ■ exp ±i\/~29 a 


* 


1 


* 





±1 


4itp-CDW 


cos 2\/2<f> p 


* 


-1 


* 


4k F 


±1 


singlet 


cos \/2(j)<j 


1 


1 


1 





=Fl 


tripleto 


sin y/2~4>cr 


-1 


-1 


-1 





Tl 


triplet±i 


exp ±i^/29 a 


* 


1 


* 





Tl 


"dimer" 


cos V24> p 


1 


1 


* 


2k F 


Tl 


"Neel" 


sin \f2<j>p 


-1 


-1 


* 


2k F 


Tl 


"doublet" 


exp ±iv20p 


* 


1 


1 





Tl 



TABLE II. Discrete symmetries of the excitation spectra 
(C: charge conjugation, V: space inversion, T: spin rever- 
sal, and k: wave number). BC= 1 (BC= —1) stands for 
(anti)periodic boundary conditions. The upper (lower) sign 
of BC denotes N/2 = odd (even) cases, where iV is the number 
of electrons. The upper 12 states are "physical" ones, which 
appear under the same BC as those of the ground state. The 
lower six states are the "artificial" ones extracted by twisting 
BC with respect to the ground state. 
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Next, we consider the symmetri es o f operators. The 
operator of the marginal field (13a) has the same 
quadratic form of the Gaussian part of the Lagrangian 
density of Eq. (^), so that it has the same symmetry 
as the ground state (C = V = T = 1, k = 0). We 
can find the symmetries of the O v \ and the U 2 opera- 
tors by considering the change of the phase fields. Since 
we restrict our attention to the Hilbert space with fixed 
electron number and total z-spin, we do not consider the 
change of the 9 U fields in the symmetry operations. At 
half-filling (k F = tt/2), it follows from Eq. (|) that the 
phase fields <j> v are transformed under particle- hole (C: 
4>r,s Ws)> space-inversion (V: R^L, x — ► x + a), 
spin-reversal transformations (T: T^i), and shift oper- 
ation (S : x — ► x + a) as 



C 
V 
T 
S 



b<r, 4> P -> -<f> P , 

t>cr, 4>p 
<f>p ■ 



Vpi 



(29a) 
(29b) 
(29c) 
(29d) 



In this case, CV = 1 is always satisfied, so that the 
independent discrete symmetries are V, T, and S. At 
quarter-filling (fcp = tt/4), the phase fields change as 



4>p -> 7T/V8 - (f) p 



(30) 



Thus the discrete symmetries of the operator O p \ for 
q = 2 arc T 5 = — 1. The relations between the operators 
and their symmetries are summarized in Table O. 

In the present numerical calculation based on the 
Lanczos algorithm, the identification is performed by pro- 
jecting the initial vector as 



|* intt ) = i(i±7>)(i±T)|i), 



(31) 



where the signs in front of the operators correspond to 
their eigenvalues, and |i) is a configuration that satisfies 
V,T\i) ^ \i). Furthermore, \i) is classified by the wave 
numbers k — 0, n. 



IV. PHASES 

In this section, we discuss the character of each phase 
that appears in the phase diagrams. In general, there 
are no long-range orders (LRO's) in ID systems due 
to strong quantum fluctuations, so that, in such cases, 
the phases are characterized by the dominant correla- 
tion functions. The correlation functions ( p"l| ) including 
the logarithmic corrections are given by integrating the 
renormalized scaling dimensions over the RG trajectory 
as 



exp 



ln(r/a) 



dl 2[ Xpi (l) 



i(0] 



(32) 



First, we consider Eq. (g|) for the spin and charge 
degrees of freedom independently. For the spin part 
which has the SU(2) symmetry, the singlet (O a i) and 
the triplet {O a 2,z) correlation functions with logarithmic 
corrections are obtained explicitly in the gapless region 



a 



R al = -ln- 3/2 (r/a), 
r 

Ra2,3 = -ln 1/2 (r/a), 



(33a) 
(33b) 



where we have used Eqs. (jl6|), @, and @. Therefore, 
the triplet correlation is more logarithmically dominant 
than the singlet one. On the other hand, when the spin 
gap opens (g* ± = -co), the singlet excitation degener- 
ates with the ground state in the thermodynamic limit, 
so that the singlet correlation becomes constant, while 
the triplet one decays exponentially. In this way, the 
triplet correlation is suppressed in the spin-gap region. 

For the charge part, at half-filling, explicit forms of 
the correlation functions including logarithmic correc- 
tions are not obtained except for the BKT or the Gaus- 
sian lines. On the BKT line, the exponent is renormal- 



ized as if* = 1, 



so that K* > 1 is always satisfied in 



the gapless region (g^j_ — 0), and then the correlation 
for the "doublet" (O p s) is dominant. In the charge-gap 
region with g% ± = oo, the "Neel" (C P 2) state degener- 
ates with the ground state, and the "dimer" (O p i) and 
the "doublet" correlations decay exponentially. On the 
other hand, for g^ ± = — oo, the "dimer" state degenerates 
with the ground state, and the "Neel" and the "doublet" 
correlations decay exponentially. 

Next, we discuss the physical states that consist of the 
charge and the spin parts. In the metallic region (g^j_ = 
0, K* > 1), the triplet superconducting (TS) correlation 
is dominant when the spin part is gapless. The operators 
for the TS phase consist of the "doublet" and the triplet 
ones, 



CtS = X] C 0S C ]+1 



exp [+iv / 26» p (x)]sin[\/2(/) (T (a;)], 



(34a) 



°T Sl =4 C kl,t> 



•j'T j'+i,T ' 

oc exp (x)) exp [+iV29 a (x)}. (34b) 

On the other hand, the singlet superconducting (SS) cor- 
relation whose operator is given by the "doublet" and the 
singlet ones 



°ss = ct t cj il 

oc exp[+iv / 26'p(j:)] cos[v / 2^ tT (x)], 



(35) 



is dominant when the spin part has a gap. 

In the insulating region, which corresponds to the fixed 
point ^ ± = +oo, the bond-spin-density-wave (BSDW) 
phaseEj characterized by 
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o 



BSDWq 



■IV V(cLrf ,c j+ i, s - + cj +1 ^T s Q s ,c jV ), 



'S3' 

s,s f 



Cbsdwz oc Sln 

(x)]sm[\/2(/) a (x)], (36a) 
Cbsdw± °c si n (x)] exp [±iV26 a (x)}, (36b) 

appears when the spin sector is gapless. On the other 
hand, the charge-density-wave (CDW) phase 



O 



CDW = (-lV^ct^g, 

s 

oc sin[v2^p(^)] cos[V2^<r(x)], 



(37) 



appears when the spin gap opens. In the CDW phase, 
both charge and spin gaps open, so that a LRO exists. 



(a) CDW 



(b) SDWz 



(c) BCDW 



(d) BSDWz 




-C__ 



-C~_ 



FIG. 3. Schematic illustration for the four charge-gapped 
states in up- and down-spin subsystems. The enclosed two 
sites in (c) and (d) stand for electron-hole dimers. The elec- 
trons polarize on sites (C = V = -1) for CDW and SDW 
states, while they polarize on bonds (C — V = +1) for BCDW 
and BSDW states. The two subsystems are synchronized 
(T = +1) for CDW and BCDW states, while they are dis- 
placed by one site (T = -1) for SDW and BSDW states. 





9z± = 


ff3X = +oo 


g*3± = -oo 


9*x = 


TS 


BSDW 


SDW 


Six = -oo 


SS 


CDW (LRO) 


BCDW (LRO) 



TABLE III. Correspondence between six possible phases at 
half-filling and fixed points in the RG analysis. For g*^ = — oo 
(g*± = 0), the spin sector is gapped (gapless). For g$ ± = ±oo 
(fllx = 0) , the charge sector is gapped (gapless) . 



In the insulating region for the opposite fixed point 
= — oo, the spin-density-wave (SDW) correlation 
characterized by 

CsDWz °C cos [y/2^ p {x)]s\n[s/2^ a {x)\, (38a) 
Osdw± oc cos [y/2<j> p (x)]exp[±\y/20 a (x)], (38b) 



is dominant when the spin part is gapless. The bond- 
charge-density-wave (BCDW) phase characterized by 

CBCDW = (-1V ^2i c js C 3+hs + c] +hs C js ), 

S 

OC COS [v / 2</>p(x)]cos[v / 2<^ (T (x)]. (39) 

appears when the spin gap opens. In the BCDW phase, 
both charge and spin gaps open, so that a LRO exists. 
The correspondence between the above six phases at half- 



filling and the fixed points are summarized in Tabic III 

At quarter-filling (q = 2), the 4fcF-charge-density wave 
(4fcp-CDW) appears when the Umklapp scattering is rel- 
evant. The operator is given by the O p \ field with q = 2, 

04CDW = ^^(x^l^X^-r^X^-^ix), 
r 

oc cos[V8cj)p(x)}. (40) 

In the rest of this section, we further clarify the differ- 
ences among the four charge-gapped states (CDW, SDW, 
BCDW, and BSDW) discussed above. For this purpose, 
we change the basis of the bosonized operators from the 
charge and spin picture to the spin up and down one by 
introducing the following new phase fields: 



bp ± 4> a , 



(41) 



where s =T, I refer to the upper and lower signs, respec- 
tively. Then, the system is interpreted as coupled spinless 
fermion systems (S = 1/2 spin chains). In this case, the 
(B)CDW and z-components of the (B)SDW operators 
are given by 

Ccdw, Csdw, oc sin(V2<£J ± sin(V2> T ), (42) 
Cbcdwi Cbsdw z oc cos 

(v^ioos^x), (43) 

where the CDW and the BCDW (the SDW Z and the 
BSDW Z ) operators refer to the upper (lower) signs in 
the right-hand sides, and sin(v2</> s ) and cos(v / 2</>s) fields 
denote Neel and dimer states of the S =1/2 sp in ch ains, 
respectively. In addition, it follows from Eqs. (|29d|) and 
(pi) that the shift operation by one site gives 



4> s — > <j> s + 7r/ \/2, 



(44) 



so that both CDW and SDW 2 states are described by two 
Neel ordered spin chains, where the two sectors (s =t, 4) 
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are synchronized in the former, while they are displaced 
by one site in the latter. On the other hand, the BCDW 
and the BSDW 2 states are given by synchronized and 
displaced dimer-ordered spin chains, respectively (see 
Fig. ||). Therefore, in the BCDW phase, the charge is 
polarized on the bonds alternatively, and the spins are 
dimerized. In the BSDW state, the charge is polarized 
on the each bond, and the spins are located on the bonds 
and remain gapless like the SDW stateEa 

The discrete symmetries discussed in Sec. Ill charac- 
terize the differences among these physical states. It 
follows from Eq. (|2^), that the spin reversal symme- 
try corresponds to whether s =t,l sectors are synchro- 
nized (T = +1) or displaced by one site (T = —1). 
Similarly, the parity and the charge conjugation distin- 
guish whether the electrons are polarized on the sites 
(C = V = -1) or on the bonds (C = V = +1). This 
interpretation for the parity is consistent with the fact 
that the 4fcp-CDW with a site LRO has the odd parity 
(V = — 1) at quarter- filling. 



V. HALF-FILLING 

Using the method explained in Sec. 0, we analyze the 
phase diagram of the EHM at half-filling. Since there 
are many instabilities, we consider the phase diagram 
separately in the charge and the spin parts by assuming 
the charge-spin separation. However, for the CDW-SDW 
transition, there is a possibility that the charge and the 
spin degrees of freedom are coupled. So we discuss two 
cases where the two degrees of freedom are separated and 
coupled. And then we identify the valid scenario from the 
comparison with the result of the strong-coupling pertur- 
bation theory. For the phase separation, which is consid- 
ered to be a first-order transition, we need an approach 
different from the one applied to the charge- and the 
spin-gap phases. We discuss the way to determine the 
phase-separation boundary and check the validity of the 
result by the strong-coupling theory. Using the above 
strategy, we also analyze the EHM with the correlated 
hopping term [Eq. (Q)]. The results are summarized in 
Fig. Ill 



A. Spin sector 

First, we determine the spin-gap ph ase b oundary fol- 

By observ- 



lowing the method explained in Sec. II A. 
ing the singlet-triplet level crossing, the phase boundary 
is found to be near the U = 2V line. Since the criti- 
cal point is almost independent of the system size (see 
Fig. ^), we can determine the phase boundary without 
any extrapolations. In order to check the consistency 
of our argument, we calculate scaling dimensions of the 
singlet and the triplet excitations from Eq. (fjj) , and con- 
firm the ratios of the logarithmic corrections. Here the 



spin-wave velocity is given by the excitation spectra for 
N a (N a ) = 1 or | rio- 1 = \m a \ = 1 as 



= lim 

L — >oc 



E(L, N, S = 1, k = 2tt/L) - E (L, N) 
2tt/L : 



(45) 



where the extrapolation is performed by the function 
v a (L) — v a (oo) + A/L 2 fkJl/L 4 , which is explained by 
x v = 4 irrelevant fieldsQEa Physically, this correction 
is related to the deviation from the linearized dispersion 
relation assumed in the TL model. Thus, the ratio of 
the logarithmic corrections can be checked as 3 : — 1 for 
the singlet and the triplet states by using the following 
relation near the critical point, 



x a i + 3a; CT2 ,3 



(46) 



Here we use the numerical data of L = 8, 10, 12 systems 
to check the scaling dimensions, and extrapolate them by 
the form A + BL~ 2 + CL~ 4 as in the same way of the 
spin velocity. As shown in Fig. ||, the extrapolated data 
become 1/2 in the gapless region. Thus, the universality 
of the transition is identified as the level- 1 SU(2) WZNW 
model. 
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0.885 
0.884 
0.883 



U/t=2 







0.01 



FIG. 4. Size dependence of the critical point for the 
spin-gap transition at U/t = 2. The system sizes are 
L = 8, 10, 12, 14. From Ref. [l|. 
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FIG 5. Extrapolated scaling dimension given by Eq. ( ^ ) 
near the spin-gap critical point at U/t = 2. The TL liquid 
theory predicts the numerical values are 1/2 in the gapless re- 
gion (V < V c ). This result shows the existence of the spin-gap 
transition. The scaling dimensions for the singlet (a; CT i) and 
the triplet (3:0-2,3) excitation spectra in the L = 12 system are 
also shown. 
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B. Charge sector 

Due to the SU(2) symmetry of the 77-pairing (pp|), a 
BKT transition takes place on the V = line for U < 
region. This phase boundary is fixed on this line for 
any strength of U < 0. In Fig. ||(a), the degeneracy 
of the "Neel" (2^2) and the "doublet" (2^3) excitation 
spectra on V — corresponds to this SU(2) symmetry. 
On the other hand, for the U > region, there appear 
two relevant level crossings as shown in Fig. ||(b). The 
one corresponds to the BKT transition due to the hid- 
den SU(2) symmetry, and the other corresponds to the 
Gaussian transition, as was explained in Sec. ||. The rest 
of the three level crossings in Fig. [3] do not correspond 
to any phase transitions, because they correspond to the 
lines ±y pi p{Q) = y p o(0) < and the line y p< p(0) = with 
y P o{0) > (Gaussian fixed line), in the RG flow diagram 
of Fig. |(b). 



2 . (a) U/t=-l 

"Neel" — >- 
''dimer" — «— 
15 ""doublet" — «- 




2.0 



1.0 



0.5 



0.0 



(b) U/t=l 




Hidden SU(2) BKT 




/ 

^ . / .-^ 








Gaussian 



-1.0 



-0.5 



0.0 
V/t 



0.5 



1.0 



FIG. 6. "Dimer" (x pl ), "Neel" (x p2 ), and "doublet" (x p3 ) ex- 
citation spectra in the charge sector vs V/t in L — 8 system 
at (a) V/t = -1 and (b) V/t = 1. In the V/t < region, 
a BKT- type transition takes place at V/t = reflecting the 
SU(2) symmetry of the 77-pairing. In the V/t > region, two 
level crossings occur due to the hidden SU(2) symmetric BKT 
and the Gaussian transitions. These three level crossings give 
the Y-shaped structure in the phase diagram. 

The hidden SU(2) BKT transition obtained by the 
"dimer" -"doublet" level crossing (x p i = x p3 ) appears 



near the U = —2V line as was predicted by the g-ology. 
The size dependence of this transition at U/t = 1 is 
shown in Fig. |^. In order to check the consistency of 
our argument, we calculate the scaling dimensions of x P i 
using Eq. uq). Here, we calculate the charge velocity us- 
ing the excitation spectra for N P (N P ) = 1 as 

E(L, N, S = 0, k = 2-k/L) - E (L, N) 



= lim 

L — >oo 



2w/L 



(47) 



Using Eqs. (|2J), we check the following relation on the 
critical line, 



Xpi + x p2 + 2x : 



P 3 



(48) 



As shown in Fig. the extrapolated data become 1/2. 
Here, the extrapolation is performed as in the same way 
of the spin-gap transition. Thus, the universality class of 
this transition is identified as a BKT type. The deviation 
from the expected value 1/2 in Fig. || stands for the effect 
of the phase separation where the TL liquid theory breaks 
down. 




-0.594 



-0.596 



-0.598 







0.01 



FIG. 7. Size dependence of the critical point of the BKT 
transition due to the hidden SU(2) symmetry at V/t = 1. 
The system sizes are L = 8, 10, 12, 14. 
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FIG. 8. Extrapolated scaling dimension given by Eq. ( |4q ) 
on the BKT critical line. The TL liquid theory predicts the 
numerical values are 1/2. This result shows the existence of 
the BKT-type transition. 

The Gaussian transition takes place along the U = 2V 
line as was predicted by the g-ology. The size dependence 
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of this transition at U/t — 3 is shown in Fig. 0. It follows 
from Eq. (27), the following relation should be satisfied 
just on the Gaussian transition line, 



Xpl + X p2 _ 1 

2 Xp3 ~ 4' 



(49) 



The result is shown in Fig. [10| The extrapolated data 
become 1 /4 from the weak- to the intermediate-coupling 
region. Thus, the transition is identified as a Gaussian 
type except for the strong-coupling region. 








l/L 1 



0.01 



FIG. 9. Size dependence of the critical point for the 
Gaussian transition at U/t = 3. The system sizes are 
L = ^yilO, 12, 14. The data agrees with the result of Cannon 



et alH that V c /t 
0.30 



1.65lg;^. From Ref. ([§ 




0.10 







1 



2 

U/t 



FIG. 10. Extrapolated scaling dimension given by Eq. ( |49| ) 
on the Gaussian critical line. The TL liquid theory predicts 
the numerical values are 1/4. Fhis result shows the existence 
of the Gaussian transition. 



C. Transition between CDW and SDW phases 

We have determined the spin-gap and Gaussian transi- 
tion lines near the U = 2V line assuming the charge-spin 
separation, however, if the g%\\ term in Eq. (^) is relevant, 
the charge-spin coupling may take place. Therefore, we 
should consider the possibility that the charge and the 
spin degrees of freedom are not separated, and that a di- 
rect transition between the CDW and the SDW phases 
takes place. To examine this possibility, we also observe 



the level crossing of excitation spectra of the CDW and 
the SDW operators (see Table which consist of both 
charge and spin components. These spectra can be ob- 
tained under conditions given in Table [n]. 



16 h++ t/A=3 

4 + + + + + CDW 

1.4 - 



CDW-SDW 



1.2 



SDW 



08 - in^**'*t:::;°' 3 
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/ singlet . 



'dimer" 



0.4 



spm-gap 
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1.3 



1.4 



1.5 

V/t 



1. 



1.7 



FIG. 11. Six excitation spectra vs V/t in L = 8 system at 
U/t = 3 near the U = 2V line. 

The level-crossing points for the three assumed tran- 
sition lines are close to the U — 2V line, but slightly 
deviate (see Fig. |ll|). The deviations from U = to 
U = oo are shown in Fig. [l2|. These three lines coin- 
cide in the weak- and the strong-coupling limits. For the 
Gaussian line, the finite-size effect is small for all regions. 
For the spin-gap phase boundary, the finite-size effect is 
small in the weak-coupling region, but it becomes large 
in the strong-coupling region. The direct CDW-SDW 
level-crossing point has large size effect for all regions. 

In order to identify the actual transition lines from 
these three lines, we use.the strong-coupling perturbation 
theory following HirschQ and van Dongen.Q The energy of 
the SDW state in the strong-coupling region of the EHM 
are analytically obtained up to the fourth order. If we 
include the correlated hopping term (|2|) in the EHM, the 
energies of the CDW and the SDW states are given by 



ffcow _ U 2(1 -0 2 t 2 
2 



L 2 (3u - 1)U 

(1 - £) 2 [(36w 2 -5v- 1)(1 - 2 - 8(3w - l)v] t 



(50) 



L 



= V - 



v{3v- l) 3 (4u- 1)U 3 
4(1 -i)H 2 In 2 



-9C(3) 



{l-v)U 

(l-£) 2 [2(1 -0 2 



(51) 



1 + v t 4 



(l-v) 3 U 3 



where v = V/U and £ = X/t. In Eq. @, we have used 
the Betjjtej-ansatz result of the 5 = 1/2 Heisenberg spin 
chain:! 



(5, • Si 



In 2, 



4 

i = -41n2 + |C(3). 



(52a) 
(52b) 
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The phase boundary between the CDW and the SDW 
phases is given by the equation^ 



(53) 



The strong-coupling theory shows a good agreement with 
the Gaussian transition in the charge part, among the 
three transition lines that we have considered. We should 
also note that the present Gaussian critical point agrees 
with Cannon et al.'s result obtained Jay the direct eval- 
uation of the CDW order parameter:li3 V c /t — 1.65^o 05 
for Ult = 3, and V c /t = 2.92 ± 0.04 for U/t = 5.5 (see 
Figs. S and O). 
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FIG. 12. Three possible transitions (Gaussian, spin-gap, 
and direct CDW-SDW transitions) along the U = 2V line of 
the EHM calculated in L = 8, 10, 12, 14 systems. The result 
of the strong-coupling expansion [Eq. ( p^ with £ = 0] agrees 
with the Gaussian transition. This means that the actual 
transitions are the Gaussian and the spin-gap transitions, and 
a BCDW state exists between them. From Ref. [l6| 

From the above results, we conclude that the ac- 
tual transition near the U = 2V line is not a direct 
CDW-SDW transition, but two independent Gaussian 
and spin-gap transitions, at least from the weak- to the 
intermediate-coupling region. In the strong-coupling re- 
gion, these two boundaries approach and appear to be 
coupled at finite U and V. Unfortunately, in the present 
analysis, we cannot determine this tricritical point, but 
it is considered to be identical to the crossover point be- 
tween the second- and the first-order transitions. This 
phenomenon is considered to be an effect of the charge- 
spin coupling term [the lo^y term in Eq. (Q)] as was dis- 
cussed by Voit in Ref. |j In this way, our analysis sug- 
gests that the crossover along the U = 2V line is closely 
related to the validity of the charge-spin separation. Our 
result also demonstrates that there is a finite region of 
a charge- and spin-gapped state between the Gaussian 
and the spin-gap transition lines. It follows from the dis- 
cussion in Sec. IV that the third phase is identified as a 
BCDW state. 

In order to clarify the above interpretation for the 
CDW-SDW transition, we analyze the EHM including 



the correlated hopping term [Eq. (0)]. This term is known 
to enlarge the BCDW phase for X/t < even in the 
U, V — > limit, without disturbing the Y- and the I- 
shaped structure of the phase diagram for the charge and 
the spin parts (see Fig. |l|). Because this term makes the 
magnitude of the backward and the Umklapp scatter- 
ing couplings different [gi± = U — 2V + 4X/ir,g3± = 
-(U-2V- 4X/tt)] conserving the SU(2)®SU(2) sym- 
metry of the Hubbard modelEj As was predicted by the 
g-ology, the BCDW phase appears from the weak- to the 
intermediate-coupling region.E3 On the other hand, in the 
strong-coupling region, the two transition lines are cou- 
pled, and the direct CDW-SDW transition takes place. 
The strong-coupling the ory also agrees with the Gaussian 
line as is shown in Fig. |l3j(a). This is the same behavior 
as in the case of the pure EHM. Therefore, the tricritical 
point in the pure EHM is considered to have the same 
property of the one that separates the CDW, the SDW, 
and the BCDW phases in the case of X/t < 0. 
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FIG. 13. The Gaussian (G) and the spin-gap (S) transitions 
in the strong-coupling region for (a) X/t = —1/4 and (b) 
X/t = 1/4 . The strong-coupling theory [Eq. (p3()] agrees 
with the Gaussian lines. Note that a BCDW phase appears 
not only in (a) but also in (b). 

For X/t > 0, the order of the Gaussian and the spin- 
gap transitions becomes vice versa, so that the BSDW 
phase appears. However, in Fig. |l3|(b), we find that the 
two transition lines cross in the intermediate-coupling re- 
gion, and a finite BCDW phase appears. This BCDW 
region becomes narrower as the system size is increased. 
In the present analysis, we cannot conclude whether 
this BCDW phase remains or vanishes in the thermo- 
dynam ic lim it. We will consider this phenomenon again 
in Sec. IvS 



D. Phase separation 

Here, we determine the phase-separation boundat 
from the numerical data of the exact diagonalization.c 
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Usually, the phase boundary-is determined by the diver- 
gence of the compressibility E3 However, for the U/t 3> 1 
region of the EHM, the phase separation can take place 
in the SDW state where the compressibility cannot be 
defined. Consequently, the method of observing the di- 
vergence of the compressibility is no longer valid. In this 
paper, we determine the phase boundary by comparing 
the energy of the ground state and that of the phase- 
separated state. In the phase-separated state at half- 
filling, the system is separated into doubly occupied sites 
and a vacuum. In this case, the energy in the thermody- 
namic limit is exactly obtained as 



Bps — — - — L. 



(54) 



Therefore, we use thep^elation E = Eps as a criterion for 
the phase separation]^ We show in Fig. [l4| the ground- 
state energy of the EHM measured from the fully phase 
separated state. The zero point gives the critical point for 
the phase separation. The finite-size effect of the phase 
boundary is sufficiently small. 
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FIG. 14. Ground-state energy of the EHM for L = 6-14 
systems at U/t = 8 measured from the fully phase separated 
state E PS = (U + 4V)L/2 [Eq. Q]. The zero point gives the 
critical point for the phase separation. 

In order to check the validity of the numerical results, 
we compare our result with the asymptotic phase bound- 
ary in the strong-coupling limit. For U/t 1 region, 
usjxug Eqs. ( pl| ) and (Q), the phase boundary is given 



E 



PS 



SDW- 



(55) 



On the other hand, for U/t <C —1 region, the system 
can be mapped onto the 5=1/2 XXZ spin chain with 
antiferromagnetic coupling J xy = At 2 /\U\ and J z = J xy + 
4 V, by using the 77-pairing operators (gOj) and the second- 
order perturbation theory. Then, the phase boundary 
is given by J xy = —J z , which corresponds to the first- 
order transition between .the XY and the ferromagnetic 
phases in the spin system.tffl If the effect of the correlated 



hopping term (0) is included, the phase boundary is given 
by 



V 



ot 2 



(56) 



The numerical result given by E = Eps well agrees with 
these asymptotic phase boundaries (see Fig. Il5j). 
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FIG. 15. Phase diagram of the ID EHM determined 
by the data of the L = 12 system at half-filling for (a) 
X/t = -1/4, (b) X/t = 0, and (c) X/t = 1/4 [CDW (SDW), 
charge (spin)-density wave; BCDW (BSDW), bond-charge 
(spin)-density wave; SS (TS), singlet (triplet) superconduct- 
ing phase; PS, phase-separated state]. The asymptotic phase 
boundaries for the PS are given by Eqs. (p5|) and (pq). 
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VI. QUARTER-FILLING 



Finally, we j 



fif the phase diagram of the EHM at 
£j by the level-crossing approach. In 



quarter-fillingcj' 
this case, the metal-insulator transition is considered as a 
BKT transition due to the higher-order Umklapp scatter- 
ing (q = 2). Then, the phase boundary should be given 
by the level crossing between the marginal and four times 
of the 4fcp-CDW spec tra x p — Ax p i or 4x p2 = x p3 , as 
was discussed in Sec. [IB 2. In the present numerical 
analysis, we use the former level crossing, because x P 3 
needs larger Hilbert space than x p o and x P 2- Based on 
this assumption, we obtain the result shown in Fig. [if]. 
Since the level-crossing point for the transition is higher 
(x p = 2) than the case of the half-filling (x p — 1/2), the 
finite-size effect from the irrelevant field (the deviation 
from the linearized dispersion relation) becomes larger. 
In this case, we need an extrapolation of the critical point 
to make the phase diagram. In the U/t — ► oo limit, 
the transition point for the charge-gap phase is given by 
V c /t = 2, because it corresponds to_th£ XY-Neel transi- 
tion in the S = 1 /2 XXZ spin cham-EH For the U/t > 1 
region, the extrapolated phase boundary flows into the 
point (U, V) — (oo, 2t) as we expected. In order to check 
the consistency in the finite-?/ region, we calculate the 
following averaged scaling dimension 



3x 



P 2 



(57) 



Except for the large- U region, the extrapolated value be- 
comes 1/2 with error less than 4%, as shown in Fig. [Tt]. 
Thus, the universality class of the transition is considered 
to be a BKT-type. 
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FIG. 16. Phase diagram of the ID EHM at quarter-filling 
for the V/t > region. The phase boundary between the TL 
liquid and the 4&;f-CDW phases is determined by the level 
crossing of x p o = Ax P 2 in L — 8, 12, 16 systems. The critical 
points in the strong-coupling limits are V c /t = 2 and U c /t = 4, 
respectively. 

On the other hand, for V/t ^ 1 region, since the 
finite-size effect is too large, it is hard to perform the 



systematic extrapolation. However, the phase boundary 
appears to flow into the exact transition point U c = At 
in the V/t — > oo limit as the system size is increased. 
Now let us review how the critical point of. the charge- 
gap phase in the V/t — > oo limit is derived. E3 The charge 
gap is defined by 



A p = E(N + 1) + E(N - 1) - 2E(N). 



(58) 



At quarter-filling, E(L/2) = 0. If one electron is add to 
this, then the energy is E(L/2 + 1) = U. Conversely, 
if one electron is removed, two free holes appear, then 
they have a kinetic energy E(L/2 — 1) ~ — At cos(ir/L). 
Therefore, the critical point for the charge-gap phase is 
given by U c = At in the thermodynamic limit. We should 
note that the critical point in the U/t — * oo and V/t — > 
oo limits are given by a completely different argument. 
Therefore, we expect that a tricritical behavior may also 
be seen along the BKT transition line. 
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FIG. 17. Scaling dimensions given by Eq. @ on the BKT 
line. The extrapolated value becomes 1/2 with error less than 
4%. 

In a similar way as the half-filling, we consider the 
effect of a charge-spin coupling operator which is de- 
rived from the q = 2 Umklapp scattering between three 
parallel spins and one antiparallel spin [the g 3 \\ term in 
Eq. (0)]. This operator may cause the synchronization 
of the BKT transition (q = 2) in the charge part and 
the spin-gap transition. In order to examine this pos- 
sibility, we determine the spin-gap phase boundary by 
the singlet-triplet level crossing in the large-U region. 
Then, the phase boundary of the spin gap appears to be 
coupled with that of the charge gap, and flows into the 
point (U,V) = (it, oo), as the system size is increased. 
Therefore, this result suggests that the crossover in the 
metal-insulator transition also takes place by the mech- 
anism similar to the case of half-filling. _Ln the large-U 
region, there are phase-separated states,E3 but that has 
not been studied in the present analysis. We should also 
consider the effect of the phase separation in the future. 
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VII. SUMMARY AND DISCUSSION 

Let us summarize the results obtained in this paper. 
We have determined the phase diagrams of the ID EHM 
at half- and quarter-filling using the level-crossing ap- 
proach, which is based on the TL liquid theory and the 
renormalization group. The metal-insulator transitions 
in half- (quarter-) filling are classified as BKT-type tran- 
sitions due to the first- (second-) order Umklapp scat- 

tization" 
and the 



tering. This fact reflects the "fracti 
discussed by the bosonization theor" 12 



generalized Lieb-Schultz-Mattis theoremx 

In the case of half-filling, for the charge sector, there 
are two BKT lines reflecting the SU(2) and hidden SU(2) 
symmetries, and one Gaussian line. These three critical 
lines meet at the multicritical point (U/t,V/t) — (0,0), 
and they form a Y-shaped structure. Note that the same 
structure is also known in the ptese diagram of the S = 
1/2 frustrated XXZ spin chainOE 2 ] For the spin sector, a 
spin-gap transition occurs due to the attractive backward 
scattering. Since the transition takes place at the origin 
of the RG flow diagram, the phase boundary has an I- 
shaped structure. Thus, the entire phase diagram is given 
by the combination of the "Y" and the "I" . 

The transition between the CDW and the SDW phases 
has been considered by assuming the following two sce- 
narios: one is the independent Gaussian and spin-gap 
transitions under the charge-spin separation, the other 
is a direct CDW-SDW transition under the charge-spin 
coupling. By checking the relations between the scal- 
ing dimensions, and by comparing the numerical re- 
sults with the strong-coupling theory, we have concluded 
that the first scenario is realized from the weak- to the 
intermediate-coupling (U ~ At) region, and the second 
scenario is realized in the strong-coupling region. In the 
former BCDW phase exists between the CDW 

and the SDW phases. Thus, the crossover of the CDW- 
SDW transitions from the second order to the first order 
turned out to be a phenomenon that reflects the validity 
of the charge-spin separation. 

To clarify the mechanism of the crossover in more de- 
tail, we have investigated the phase diagram of the EHM 
including the correlated hopping term (^|). In this case, 
there appear BCDW and BSDW phases, depending on 
the order of the Gaussian and the spin-gap transitions. 
Note that the direct CDW-SDW transition also takes 
place in the strong-coupling region. 

Therefore, we can understand the crossover of the 
CDW-SDW transition in the EHM as a kind of these 
generalized cases. The reason why the mechanism of 
the CDW-SDW transition has been left ambiguous for 
long times is that the analytical solutions both in weak- 
and strong-coupling limits give the same phase boundary 
U = 2V, and the numerical analysis for the intermediate- 
coupling region does not have enough precision to distin- 
guish the Gaussian and the spin-gap transition lines. Re- 
cently, the similar mechanism of transition was reported 



in studies of the Hubbard chain with periodic poten- 
tial, whic h ha s the transition between Mott and band 
insulatorsMa 

At quarter-filling, the phase boundary of the metal- 
insulator transition has also been determined by assum- 
ing that it is the BKT-type transition of the higher-order 
Umklapp scattering. The obtained phase boundary is 
consistent with the known exact results in the U — > oo 
and the V — > oo limits. In the present parameter space, 
there are neither a BKT transition with the hidden sym- 
metry, nor a Gaussian transition. Although the finite-size 
effect in V/t ^S> 1 region is large, the phase boundaries of 
the charge- and the spin-gap phases appear to be coupled 
as the system size is increased. Therefore, there may be 
a crossover from the BKT to the first-order transitions 
in this region. 

The rest of the section is devoted to discussions. We 
have clarified the mechanism of the CDW-SDW transi- 
tion of the EHM at half-filling. However, we have not re- 
vealed the reason why the intermediate state in Figs. [l2| 
and |J|(b) is not a BSDW but a BCDW. We can interpret 
the appearance of the BCDW phase considering roles of 
the i?3j| term in Eq. (^). When the charge gap opens in 
the U > 2V region (g^x — ~°°)7 t ne phase field is locked 
as (j) p — 2nn / '\/8 with n being an integer, which min- 
imize the classical potential energy associated with the 
<?3_l term. In this case, the charge part of the g%\\ term is 
also locked, so that gix — U — 2V is reduced as 



9t± -> 9i± +5311 (cos 



(59) 



where g$u = — 2V. Thus, if we estimate the spin- 
gap pha se boundary in the weak-coupling region using 
Eq. @, it shifts toward the U > 2V side of the U = 2V 
line. Similarly, we can estimate the shift of the Gaus- 
sian line. In the U < 2V region, the spin gap opens 
(g*j_ = —00), then g%± = —U + 2V shifts as 



93± -> 93± + .9311 (cos[V8</v]), 



(60) 



where the phase field is locked as <\> a — 2n7r/v / 8- There- 
fore, the Gaussian line shifts toward the opposite side of 
the spin-gap phase boundary. Consequently, the BCDW 
phase appears between the CDW and the SDW phases. 
Thus, it turns out that the <? 3 || term enhances the BCDW 
phase when it is irrelevant, and it couples the Gaussian 
line and the spin-gap phase boundary when it is rele- 
vant. The deviations from the U = 2V line may be an- 
alyzed quantitatively in the weak-coupling region by the 
renormalization group analysis including g^n. The above 
explanation may also be applicable to the phenomenon 
that a BCDW phase appears for X/t > in the strong- 
coupling region [see Fig. |l3|(b)]. 

We should consider the reason why the size depen- 
dence of the Gaussian transition line is smaller than that 
of the spin-gap phase boundary, in the strong-coupling 
region of Figs, [l^ and [l3| The reason is considered to 
be the difference between the behavior of the charge and 
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the spin gaps. As was discussed in Sec. & the spin gap 
opens exponentially slow (jl8|) , while the valley of the 
charge gap near the Gaussian transition becomes steeper 
as the strength of the interaction is increased ( ^8|) . In the 
present analysis based on the TL liquid theory, we have 
assumed that both charge and spin parts are gapless, so 
that the Gaussian (spin-gap) transition line may be af- 
fected by the spin (charge) gap in the strong-coupling 
region. In the present case, magnitude and variation of 
the charge gap are considered to be much larger than 
those of the spin gap near the U = 2V line in the strong- 
coupling region, so that the Gaussian line has less size 
dependence than the spin-gap phase boundary. 

In the present paper, we have not determined the 
tricritical poinL. || Th is.|-problem has been discussed by 
many authors. BOO'tLSta Recently, the tricritical point 
is explo red b y the density-matrix-renormalization-group 
method,t§E3 but we may also determine the tricritical 
point by comparing the numerical result of the Gaussian 
transition line and the strong-coupling perturbative ex- 
pansions in the higher order. It may also be worth study- 
ing the level crossing for the tricritical point by consider- 
ing the logarithmic corrections to the excitation spectra 
which stem from the term. In addition to the tri- 
critical point between the CDW and the SDW phases, 
there is another tricritical point that separates the TS, 
the SDW, and the PS states. In Fig. [l5|, the phase bound- 
aries between the SDW and the PS states obtained by the 
strong-coupling theory seem to be bent near the tricriti- 
cal point. This tricritical point may also be identified by 
the further strong-coupling calculation. 

Finally, we discuss the effect of site-off-diagonal inter- 
actions. In this paper, we have considered the correlated 
hopping term given by Eq. however the generalized 
form of this term is given byEHI 

Hx = X ^(c\ s c t+lyS + H.c.)(n;,_ s + m + i- s ), (61a) 

is 

Hx> = X' ^T(cl s Ci + i >a + R.c.)m,- s n i+ i _ s . (61b) 

is 

In the present paper, we have set X = —X'/2 to keep the 
particle-hole symmetry and the SU(2)<g)SU(2) symmetry 
of the Hubbard model. If this relation is not chosen, 
these symmetries are lost, so that the V = line is no 
longer the phase boundary of the metal-insulator transi- 
tion at half-filling. Besides, an additional UmklanpLterm 
sin \/8(j)p appears in the effective Hamiltonian.QEj The 
analysis for this situation (X ^ —X'/2) is the subject of 
future research. 

There are other types of site-off-diagonal interactions. 
Foj^example, the bond-bond interaction term is given 
byN 

H w = ^^(4q +li8 +H.c.)(cU +1 , 8( +H.c). (62) 

iss f 

This term also makes the difference in the magnitude 
of the g-parameters for the backward and the Umklapp 



scattering in the weak-coupling limit \g±± — U — 2V + 
8 W, g 3J _ = - (U - 2V - 8W)] (Ref . |) conserving the sym- 
metries of the Hubbard model. Therefore, we expect that 
this term play a role similar to the correlated hopping 
terms [Eq. (g)]. On the other hand, since this term can 
be rewritten as the exchange of the spins and of the pseu- 
dospins ([?0]), it may_affect the first-order transition to a 
ferromagnetic stateEa or the phase separation. Besides, 
in a parameter region of the EHM including this term, 
the BSDW state is shown to be the exact ground state.E^I 
The analysis of the effect of this term hv the level-crossing 
approach will be reported elsewhereE3 
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APPENDIX A: WEAK-COUPLING LIMIT 

In the weak-coupling limit, the parameters of the 
sine-Gordon model (^) can be identified in terms of 
the bare coupling constants of the original model, 
and consequen±Wr-the phase boundaries are obtained 
j'BB'aEJ This approach is often called the g- 



analyticallym 
ology. 

The parameters of Eq. 
parameters defined in Rcfs. 



, are given by the 
]3| and [| as follows 



'9u_ 
.2tt, 

9411 ± 54-L 



2tt 



K v = 
9u = 



2-ku v + g v 
2-nu v - g v ' 

.9i|| - 3211 Tff2_L, 



(Al) 



where the upper (lower) sign corresponds to v = p (v = 
a), and vp — 2tsin(fcp) is the Fermi velocity. For the 
EHM including the X [Eq. (§)] and W [Eq. @] terms at 
half- filling, the g-parameters can be identified as follows: 



9i± = g a = U ~2V + 6g, 

g»± = -(U-2V-8g), 
g p = -(U + QV-Sg), 
Sg = 4X/n + 8W. 



(A2) 



This calculation can be performed straightforwardly ex- 
cept for the X term. For the X term which contains a 
three-body term, the operator-product-expaasion tech- 
nique is needed to identify the g-parametersH3 

The instabilities for the charge and the spin gaps are 
discussed based on the renormalization group (RG) anal- 
ysis as explained in Sec. y. Since K u is approximated as 
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K u w 1 + g v /2-Kv v , the parameter in the RG flow diagram 
(Fig. |) is given by y 0v (l) = g v /irv v . Then, the spin-gap 
opens when gi± = g a < 0, and the phase boundary is 



U = 2V-Sg. 



(A3) 



On the other hand, the charge gap opens when g$± > 
\g p \. The condition 93^ = — g p < is the BKT-type 
transition due to the SU(2) symmetry in the charge sec- 
tor. Then, the phase boundary is obtained as 



V = Q, U < 



(A4) 



The condition g%± — g p > is the BKT-type transition 
due to the hidden SU(2) symmetry in the charge sector. 
Then, we obtain the phase boundary as 



U = -2V + (5 5 , U > Sg. 



(A5) 



The Gaussian-type transition takes place at g%± = 0, for 
g p < 0. Thus, we obtain the Gaussian line for the charge 
sector as 



U = 2V + Sg, V < 0. 



(A6) 



Thus, we obtain the phase diagrams in the weak-coupling 
limit as is shown in Fig. [I] which have the Y- and the I- 
shaped structures in the charge and the spin degrees of 
freedom, respectively. In the present g-ology analysis, 
the parameters X and W appear only through 5g, so 
that the X and W terms play similar roles in respect of 
the enhancement of the BCDW (Sg < 0) or the BSDW 
(Sg > 0) phases. The correspondence between the g- 
ology and the level-crossing approach is summarized in 
Table |. 

In this paper, we have taken U rjS into account in the 
identification of g\± and g^± by following Ref. The 
Clifford algebra for U r , s @ can be expressed by tensor 
products of the Pauli matrices, which are chosen so the 
effective Hamiltonian (^|) is to be diagonal in the space 
of U r , s - For example, the following representation is pos- 
sible: 



T 



T 



1 



(A7) 



These operators contribute to the gn_ and the g%±^ terms 
as ±1 (g) t z , and one eigenvalue of the matrix is chosen. 
Consequent ly, t he signs of <?u_ and g%± become opposite 
as in Eqs. (A2). When we consider the bosonization of 
a physical operator, it should be diagonalized simultane- 
ously with the Hamiltonian. Therefore, if t he physical 
operator cannot be diagonalized by Eq. (A7), we should 
choose other representations for U r s . In the derivation 
of the physical operators in Sec. |v|, we need one more 
representation such as 



E/rt = t z 
C/ LT = r y 



1 r- 



u u = 1 



(A8) 



In this case, the contribution to the gi± and the g3± 
terms is =Fl ® t z . If the contribution of U TtS is neglected, 
the sign of the <?3j_ term is reversed, and the roles of O p \ 
and P 2 are interchanged. 
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